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Abstract 

We characterize completely the Gneiting class [6] of space-time covariance functions and give more relaxed 
conditions on the involved functions. We then show necessary conditions for the construction of compactly 
supported functions of the Gneiting type. These conditions are very general since they do not depend on the 
Euclidean norm. Finally, we discuss a general class of positive definite functions, used for multivariate 
Gaussian random fields. For this class, we show necessary criteria for its generator to be compactly support 
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1 Introduction 

Recent literature persistently emphasizes the use of approximation methods and new methodologies for dealing 
with massive spatial data set. When dealing with spatial data, calculation of the inverse of the covariance matrix 
becomes a crucial problem. For instance, the inverse is needed for best linear unbiased prediction (alias kriging), 
and is repeatedly calculated in the maximum likelihood estimation or the Bayesian inferences. Thus, large spatial 
sample sizes traduce into big challenges from the computational point of view. 

A natural idea that made proselytes in the last year is to make the covariances exactly zero after certain distance 
so that the resulting matrix has a high proportion of zero entries and is therefore a sparse matrix. Operations on 
sparse matrices take up less computer memories and run faster. However, this should be done in a way to preserve 
positive definiteness of the resulting covariance matrix. The idea goes under the name of covariance tapering, 
by meaning that the true covariance is multiplied pointwise with a compactly supported and radial correlation 
function. This operation is technically justified by the fact that the Schur product preserves positive definiteness. 

The effects of tapering in terms of estimation and interpolation have been recently inspected by [4] , where general 
conditions are given in order to ensure that tapering does not affect the efficiency of the maximum likelihood 
estimator. For spatial interpolation, [5] show that under some regularity conditions, tapering procedures yield 
asymptotically optimal prediction. In order to assess these properties, the asymptotic framework adopted by the 
authors is of the infill type, and the tool allowing to evaluate the performances of tapering is the equivalence of 
Gaussian measures, for which a comprehensive theory can be found in the seminal work by Yadrenko |13j . 

These points fix very briefly the state of the art and we refer the reader to [3] for an excellent survey on the 
topic. 

Although tapering has been well understood in the spatial framework, there is nothing done, to the knowledge 
of the authors, for the spatio-temporal case. In particular, the use of tapering is at least questionable in space-time, 
since the same type of asymptotics does not apply and thus it is not easy to evaluate its performances. 

But a deeper look at this problem also highlights the non existence, in the literature, of space-time covariance 
functions that are compactly supported over space, time or both. These facts motivate the research documented 
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in this manuscript. 

We deal with challenges related to space-time covariance functions. If spatial data set can be massive, one can 
imagine how the dimensionality problem affects space-time estimation and interpolation. This problem may be 
faced on the base of two perspectives that can be illustrated through the celebrated T. Gneiting class of covariance 
functions [6]: for (x,t) £ R d+l , the function 



is positive definite, for tp completely monotone on the positive real line and h a Bernstein function. For 1 = 1, the 
function above is a stationary and nonseparable space-time covariance. This function has been persistently used by 
the literature and a Google scholar search highlights that currently there are over 90 papers where this covariance 
has been used for applications to space-time data. 

If there are many observations over space, time or both, then the use of this function would be questionable 
for the computational reasons exposed above. A more intriguing perspective is to consider a function of the 
Gneiting type, but replacing the generator ip in equation (J]) with a compactly supported function, and inspecting 
the conditions ensuring that permissibility is preserved on some d-dimensional Euclidean space. The results are 
illustrated in the following sections. 

An auxiliary result of independent interest is also given: we characterize completely the Gneiting class and give 
more general conditions for its permissibility. 

The ratio mentis of this paper leads then to consider a general class of covariances, originally proposed in 
Porcu et al. [10] and more recently in pQ. Both groups of authors show that is class of covariances can be very 
versatile since it can be used for two-fold purposes: on the one hand, it can be effectively used to deal with zonally 
anisotropic structures, on the other hand it can be adapted to represent the covariance mapping associated to a 
multivariate random field, which is highly in demand since there are very few models with these characteristics [7]. 

As a conclusion to the preludium, the plan of the paper is the following: in Section 2 we present basic facts 
about positive and negative definite functions. Section 3 characterizes completely the Gneiting class, for which 
only sufficient conditions were known until now. In Section 4 we present necessary conditions for compactly 
supported covariances of the Gneiting type. Similar results are obtained in Section 5 for the multivariate class of 
cross-covariances proposed in |10| . 

2 Preliminaries 

This section is largely expository and contains basic facts and information needed for a self-contained exposition. 
We shall enunciate the concepts of positive and negative definiteness, as well as the material related to them, 
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working with linear spaces and subspaces. The space-time notation will be used only when necessary for a clearer 
exposition of results. 

For E a real linear space, we denote by FD(-E') the set of all linear finite-dimensional subspaces of E. If 
dim £ = n£H and e\ , . . . , e„ are basis in E, then 

feC(E) ^ f{x 1 e 1 + ...+x n e n )&C{R n ) 
and 

feL(E) ^ /( Iiei + ... + Vn )GL(n 

Also, we call Cq(E) the set of all function / e C(E) such that / has compact support. If dimE = oo, then 
feC(E) ^ f e C(E ) VEq € F~D(E). 

A complex- valued function / : E — > C is said to be positive definite on E (denoted hereafter / € &{E)) if for 
any finite collection of points {£i}?=i € E the matrix (/(& — £j))I\'=i ' ls positive definite, i.e. 

n 

for all ai, a%, . . . , a n e C : ^ - £j)aj > 0. 

It is well known that the family of positive definite functions is a convex cone which is closed under addition, 
products, pointwise convergence and scale mixtures. Briefly, we have the following properties. 
Let /, fi e $(E), i e N. Then: 

1. \f(x)\ < /(0), 7F^) = fix), \f(x) - f(h)\ 2 < 2/(0)Re(/(0) - f(x - h)), x, h e E; 

2. Xih + A2/2 with Ai > 0, /, Re/, hh G 

3. if, for all x € E, the finite limit lim /„(cc) =: exists, then g € &(E); 

n — >oc 

4. for any linear operator A : E\ — > E the function / o A belongs to $>(E\); in particular, / G for any 

linear subspace i?i from E. 

Let £ = R". The celebrated Bochner's theorem establishes a one to one correspondence between continuous 
positive definite functions and the Fourier transform of a positive and bounded measure, i.e. f(x) = F n (fi(u))(x). 
If /i is absolutely continuous with respect to the Lebesgue measure, than d[i(u) = f(u)du, for / nonnegative. This 
can be rephrased in the following way: if / € C(R") n L(R n ), then / e $(R n ) if and only if 



f(u) = F- 1 (f)(u):= e i ^f(x)dx>0, u e I", 
it" 

for (•, •) the usual dot product. The function / is called spectral density or Fourier pair associated to /. 



3 



If / is a radially symmetric and continuous function depending on the squared Euclidean norm || • Hf, i.e. 
f(x) = ydlxlH), p £ C[o,+oo) ; then the Fourier transform above simplifies to the Bessel integral (if in addition 

/ G L(R")) 

/■+00 

9n(s)~ i P (u 2 )u n - 1 j%-i(su)du, (2) 

Jo 

where j'a(u) := J ^"' > , with Ja a Bessel function of the first kind. Thus / G $(R"), for some n G N and for / 
radially symmetric, if and only if g n (u) > Vu > 0. 

A function / :]0, oo[— > R is called completely monotone, if it is arbitrarily often differentiable and 

(-l)'7 (n) (x) > for x > 0, n = 0, 1, . . . . 

By Bernstein's theorem the set M( 0jOO ) of completely monotone functions coincides with that of Laplace transforms 
of positive measures \x on [0,oo[, i.e. 

f(x)=Ln(x)= I e-^dfiit), 

J[0,oo[ 

where we only require that e~ xt is ^-integrable for any x > 0. A/(o,oo) is a convex cone which is closed under 
addition, multiplication and pointwise convergence. 

The connection with the function g n {-) gives the celebrated Schoenberg (1939) theorem by which a radial 
function f(x) = ^(H^Hl), f & C[o,+oo)> belongs to $(R") for all n G N if and only if ip is completely monotone on 
the positive real line, and in this case the Bessel integral in equation ^ reduces to a Gaussian mixture. Finally, a 
Bernstein function is a positive function that is infinitely often differentiable and whose first derivative is completely 
monotone. For a more detailed exposition on these facts the reader is referred to [TT] . 

In this paper we shall be also dealing with functions depending not on the Euclidean norm but on some 
homogeneous continuous function p : E — > R such that p{tx) = \t\p(x) Vi G R, x G E and p(x) > 0, x ^ 0. If 
ip G C^+oo) and J +O ° |v5(i 2 )|i n_1 dt < +oo, then we have that tp o p 2 G <I>(R n ) if and only if the function 

i"3" H G n (v) := / V (p 2 (y))e^ dy (3) 

is nonnegative for all v G R n . If p is the Euclidean norm, then the functions G„(-) and g n {-) are related by the well 
know equality G n {v) = (27r) ^ (| |w| I2) - 

Finally, a complex- valued function h : E — > C is called (conditionally) negative definite on E (denoted h G N(E) 
hereafter) if the inequality 

n 

^ c k c 3 h(x k -Xj) < 
is satisfied for any finite systems of complex numbers ci, C2, c„, X)fc=i c ^ = ®> an d points 
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Let {Z(£),£ G R™} be a continuous weakly stationary and Gaussian random field (RF for short). The associated 
covariance function / : R™ — > R is positive definite. This can be rephrased by saying that positive definiteness of a 
candidate continuous function / : R™ — > R is sufficient condition for the existence of a continuous weakly stationary 
and Gaussian RF having /(•) as covariance function. 

If, additionally, /(•) is radially symmetric, the associated Gaussian RF is called isotropic. Isotropy and sta- 
tionarity are independent assumptions but throughout the paper we shall assume both in order to keep things 
simple. 

To complete the picture, the variance of the increments of an intrinsically stationary Gaussian RF is called 
variogram. For two points of R™, say i = 1,2, we have that Var(Z(^ 2 ) — := 7(^2 — The mapping 

7(-) : R" — > R is conditionally negative definite. The additional property of isotropy is then analogously defined as 
before. 

3 Complete Characterization of the Gneiting class 

Lemma 1. 



Proof. 

i. The necessity is obvious. As for the sufficiency, for n G N and xi, x n in E, we have that x%, ...,x n G Eo - the 

linear span of these elements. Obviously dimi^o < n. 

ii. Again, the necessity is obvious. For the sufficiency, let ei, . . . , e n be basis in E. Then we take g(xiei + . . . + 

x n e n ) = (1 — e|a;i|) + • . . . • (1 — e|x„|) + and e { 0. The proof is completed. 



i. / e HE) 



/e$(E )V£ eFD(£). 



ii. If dim E 



neN then f e $(£?) 



fg g V 5 G $(E) n C (E). 



□ 



Lemma 2. Let the next conditions be satisfied: 



1. h,be C(E) and h(t) > Vt G E. 



2. ip G Cm+oo) and for the some m G N : J Q °° \ip(u 2 )\u m 1 du < +00. 



3. p G C(R m ), p(fcc) = |t|p(a;) Vi G R, a; G R m and p{x) > 0, x ^ 0. 
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with G m (-) defined in equation (B). 

Proof. Observe that tp (p 2 {x)) G L(W n ). We have that 

K{x,t) G $(K m xE) K(x,t) G <J>(M m x E Q ) V£ G FD(E) 

K{x,t)g(t) G $(R m x S ) V£ G FD(£),V<? G ®(E ) n C (£o) 
^ // K{x,t)g{t)e l{x ^e^' u Uxdt > 

V£ G FT)(E),Vg G $(£ ) n C (£ ) , Vw G K m ,u G £ - 
As for the last integral, a change of variables of the type x = y/h(t)y yields that the last inequality is equivalent to 

g (t)b(t)(h(t))^G m (y^Mtjv)e l ^ u) dt > , Vu G l m ,iie E , 
which holds if, and only if Vg G $(E ) n C (E ), Vv G M m , we have 

g(t)b{t)(h(t))^G m (y/h^jv) G *(£b) V£ G FD(£), 
b{t){h{t))^ G m {^/h(t)v) G $(E) Vu G M m . 

The proof is completed. □ 

The following result gives a complete characterization of the Gneiting class, with the additional feature that 
only negative definiteness of the function h is required, whilst Gneiting's assumptions are much more restrictive as 
it is required that h! is completely monotone on the positive real line. Furthermore, we give a simple proof of this 
result and we defer it to next section for the reasons that will become apparent throughout the paper. 

Theorem 3. Let h G C(E), h(t) > Vi G E. Let d G N. The following statements are equivalent: 

1. K(x, t) := (h(t))-i<p G x£)V^£ C [0l+ oc) f] M (0 , + oc) ■ 

2. e - Xh{t ) G <P(E) VA > 0. 

Let us consider examples of functions h for which the statement 2 in Theorem [3] holds. 

Example 1 Let h(t) = ||£||° + c, c > 0, < p < +oo, a > 0, t = (t u ...,t n ) G W 1 , where \\t\\P = YT k= i M p , 
< p < oo, and \\t\loc = sup 1 < fc <„ |i fe |. Then 

e -xh(t) G $ ( R «) VA > ^ e-^p G $(R n ) ^> < a < a(Z") , 

2 ifn=l,0<p<oo 

p if 7i > 2, < p < 2; 

1 if n = 2, 2 < p < oo 

if n > 3, 2 < p < oo. 



where 



, p if n > 2, < p < 2; . . 

1 p J 1 1 if n = 2, 2 <p < oo; W 



For < p < 2, we get Schoenberg's result. The other two cases have been investigated by Koldobsky [8] in 
1991 and Zastavnyi [H CGH [16] in 1991 (2 < p < oo, n > 2). Finally, Misiewiez [9] gave the last result in 
1989 {p = oo, n > 3). 

Example 2 If p(t) is a norm on R 2 , then e~ p "^ G $(M 2 ) for all < a < 1. This is a well-know fact (see, for 
example, [17]). Therefore e- A ' l W G $(R 2 ) VA > 0, where /i(t) = p a {t) + c, < a < 1, c> 0. 

Example 3 Let ^>(s) 6 i Vs > 0. Then, it is well known that 

e"^' G A/(o,+oo) VA > <=> 0' G M (0 , +oo) . 

Gneiting [6] proves the following: if ip G C[ ,+oo), ip(s) > Vs > 0, and V' G M( , +oo ), then e ~ A,l(t) G $(R") 
for all A > 0, n G N, where h(t) := "0(1 1*[ li) an d, hence, 

Example 4 By celebrated Schoenberg's Theorem [12], if t) = Vt G _E, then 

G N{E) e~ xh(t) G VA > 0. 

Example 5 Let g G $(E), g(-t) = g{t) for alU G E and h(t) := g(0) - g(t) + c, c> 0. Then h(t) > Vi G £, 
h G N{E) and, hence, G 3>(E) for all A > 0. 

4 Necessary conditions for compactly supported functions of the Gneit- 
ing type 

From now on let us write § d_1 := {x G R d : ||x||2 = 1} for the sphere of M. d . 
Theorem 4. Let the next conditions be satisfied: 

1) he C{E), h(t) >0VteE and h(t) £ h(0) on E. 

2) <pe q ,+oo), <p(o) >o. 
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3) FordeN,pe C(R d ), p(tx) = \t\p(x) Vi G R,x G R d and p(x) > 0, x ^ 0. 

4) K(x, t) -.= {h(t))-iip e m d x £0- 



T/ien: 

1. {h{t))-i G <f>(E) andip(p 2 (x)) G $(M d ). 

2. if i/iere exists an £ Nf][l,d] such that / + °° \f{u 2 )\u n ^ 1 du < +oo, iften Vra = l,...,n and u G R m ifte 

function s i— > f m ,v(s) ■= s m ~ G m {sv) , with G m () as defined in is decreasing on (0,+oo). Furthermore, 
fm,v{+oo) = /or u 7^ 0. 

3. 7/ J +o ° |<^(M 2 )|M d_1 (iw < +oo, £/ien Gd(0) > 0. //, m addition, Gd is real- analytic, then Vu G R d , d ^ i/ie 

function s \— » fd, v {s) ■= Gd(sv) is strictly decreasing on [0, +oo) and Gd{v) > Vu G R d . 

4. 7//+°° | V (u 2 )K +1 dn < +oo, ftenaiH := / Rd ^(p 2 (y))(y, «) 2 dn > to £ S^ 1 and A := J Rd </>(p 2 (y))|M|! dy > 

0. Furthermore, a x (y) = on § d " 1 ^> ft = 0. If, in addition, ft > 0, then e ~ xh{t) G $(E) VA > 0. 

5. If \ip{u 2 )\e EU du < +oo /or some e > f/or example, when ip has compact support), then Bp G N : 

e -\h (t) g <j)(^) \/A > 0. In practice, for p it is possible to take one of the following numbers: 

p{v) := min { k e N : a k (v) = [ ^{p 2 (y))(y, v) 2k dy ± ol , v G S^ 1 , 



g :=min fcGN: ft = / V {p\y))\ \y\ \f dy ^ 
T/ie function p( ) is bounded on W 1 " 1 and q = min p(v). 

Proof. The statement 1 is obvious. 

Let us prove the statement 2. By Lemma El we have 

F m , v (t) := (/i(0)^G,„( V / Mi)«) G , Vm = L~n , w G 



Hence, F m ,„(0) = (/i(0))^G m (y / X(0)tO > and \F m , v {t)\ < F m)t) (0), * G £7. Therefore G m (u) > 0, v G R m , and 



(sh(t)) — G m {^/h(J)sv) < (sh(0)) — G m (^/h{0)sv) , Vm = l,n , v G R m , s > , t £ E . 
The latter inequality is equivalent to 



• «j < f m ,v 0) , Vm = l,n , v G R ro , s > , t G £. 
Since (7i(t))~3 G then ft(t) > /i(0), t £ E. Since ^ /i(0) on £7, then there exists a point t G E 



such that g := J -j^y- > 1. By the intermediate values Theorem Va G [l,gj 3£ G -E : y 7^ = a - Therefore, 



fm,v(as) < f m ,v(s) for all s > and a G [l,q\. Hence, / m , t ,(a 2 s) < f m , v {as) < f m ,v(s) for all s > and a G [l,q]. 
Thus, f m ,v(& p s) < f m ,v{s) for all s > 0, a G [1, g] and peN. This implies that the function f m ,v{s) decreases in 
s G (0, +oo). By the Riemann-Lebesgue Theorem, it follows that G m (v) — > as ||«||2 — > +oo. Hence / m , t .(+oo) = 
for i> ^ 0. The statement 2 is proved. 

Let us prove the statement 3. i. From statement 2 it follows that for all v G M. d , n ^ 0, the function Gd(sv) 
decreases in s G [0, +oo) and, hence, < Gd{v) < G<j(0). Therefore, Gd{0) > (otherwise Gd{v) = on R d =4> 
ip(p 2 (y)) = on R d , that contradicts the condition ip(Q) > 0). 

ii. If, in addition, Gd is real-analytic, then Vw G M d , !)/0, the function Gd(sv) strictly decreases on [0, +oo). 
This can be proved by contraddiction. Let us assume that, for some v G R d and vo ^ 0, the function Gd(svo) 
is constant on some interval (a,/3) C (0, +oo), a < (3. This would imply that Gd it is constant on [0, +oo) and 
Gd(0) = lim s ^ +00 Gd(svo) = 0, which contradicts i. Thus, Vv G K d , u ^ 0, the function Gd(sv) strictly decreases 
on [0, +00) and, hence, Gd{v) > lim s ^ +00 Gd{sv) = 0. The statement 3 is proved. 

Let us prove statement 4. Let v G S d_1 and fd,v(s) := Gd(sv). From statements 2 and 3, it follows that the 
function fd,v( s ) decreases on [0, +00) and that fd,v(0) > 0. Obviously, fd,v( s ) G C 2 (IR) and 



where fd. v {0) = —ai(v). Note that f'J v (0) < 0, otherwise the function fd,v(s) strongly increases on [0, c] for some 
c > 0, which contradicts statement 2. Thus, > for all v G S d_1 . For p > 0, the next integral is constant on 



Js-*- 1 

where da, if n > 2, is the surface measure on S d_1 and d<r(v) — S(v — 1) + <5(v + 1), if d = 1 (here - the Dirac 
measure with mass 1 concentrated in the point v = 0). Therefore, 



/«L(°) 



s 2 + o(s 2 ) , s -s- 0, 



2 





(5) 



Hence 




and ai(v) = on S' 



0. 



Let, in addition, /3i > 0. Then f^ v (0) 




cti(vo) < for some vq G § d 1 and 

_ \ Tl 

— ) = (1 + <?«(*))" G . Vn G N , 7« > 0. 



(6) 



Take 
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Obviously, j n — > +0 and 

9n{t) = 7 77^ 7T f TTyT ' 17" V h(t)) = h(t) , 71 -» OO . 

Therefore, Vn(*) — * e - ^*) and, hence, e _A/l W g $(E) for all A > 0. The statement 4 is proved. 
Let us prove the statement 5. In this case Gd is real-analytic and 

/^ fe) (0) - (-l)*a*(«) , /£* _1) (0) = , / a^)^) = c rf , 2fc /3 fe , fc G N. (7) 
Therefore, Vv 6 S d_1 3p G N so that 

f (2p) (0) 

fdM = U AO) + -j^r * 2p + o(s 2 n ,s^o, 

where f d 2 ^(0) ^ 0, otherwise the function fd,v(0) = fd,v(s) = /d,u(+oo) = which contradicts the inequality 
Gd(0) > (see statement 3). Hence, f^\o) < 0, otherwise the function fd,v{s) strongly increases on [0,c] for 
some c > 0, which contradicts statement 2. Thus the function p{v), v £ S d ~ 1 , defines correctly. 
Let v £ S'' -1 and p = Take function |(6]), where «o = v 

7 -r^r'»J >0,A>0 ' 

Then g n (t) £ ■ n -> oo. Therefore t/>n(*) -> e"^*) and, hence, e" A ' lP (*) £ $(E) for all A > 0. 

If ak(vo) ^ for some vq E §> d ~ 2 , k £ N, then afe(w) 7^ in some neighborhood of a point vo and, hence, 
p(v) < p(vo) in this neighborhood. Thus the function p(v) is locally bounded on compact S d_1 and, hence, p(v) is 
bounded on § d_1 . 

Let to = min p(v) = p(vq) for some Vq eS^ 1 . Then a m (v ) ^ and for all v £ S^ 1 equality 

f (2m) (0) 

fd,v(») = Uv(Q) + -^j^r s 2m + o(s 2m ) ,s^o 

(ZTOj! 

holds. Obviously {-l) k a k (v) = f d ^ ] (0) = 0, for all 1 < k < m (if m > 2), and (-l) m a„» = f£™\o) < 
(otherwise the function fd,v(s) strongly increases on [0,c] for some c > that contradicts a statement 2). From 
© follows that P k = for all 1 < k < m (if to > 2) and (-l) m /3 m < 0. Therefore q = m. 

The Theorem 2] is proved. □ 

Proof of Theorem H If h{t) = h(0) > on E, then the implication 1) => 2) is obvious. If h(t) /j(0) on 
then this implication follows from statement 4 of Theorem |4] for tp(s) = e~ s £ C^.+oo) Pl^o.+oo)- 
The reverse implication 2) => 1) follows from Lemma [2] for y?(s) = e _s , equality 

e -&\\v\\l e *lw) dy = {2 %a )i e -flHll , w G R d , > 0, 



and Bernstein-Widder's Theorem. The proof is completed. □ 
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Next Theorem [5] is an addition to Theorem |4] for the case p{x) = \\x\\i. 

Theorem 5. Let the next conditions be satisfied: 

1) he C(E), h{t) >0VteE and h(t) ^ h(0) on E. 

2) ^eq ,+oo), v(o) >o. 

3) K(x,t) := (h(t))-i(p (^f ) G $(M d x £). 

// J +o ° |<p(u 2 )|u m_1 du < +oo /or some natural m G [l,o!] and g m is real- analytic, then the function f m (s) := 
s m ~ g m (s) strictly decreases on (0, +oo) and g m (s) > for all s > 0. 

Proof. From Theorem [4] it follows that f m decreases on (0, +oo) and f m (s) > f m (+oo) — for s > 0. Since f m is 
real-analytic on (0, +oo), then function f m (s) strictly decreases on (0, +oo). Otherwise the function f m is constant 
on some interval (a, ff) C (0, +oo), a < j3, and, hence, it is constant on (0, +oo) and / m (s) = f rn (+oo) = 0, s > 0. 
Therefore, G m (v) = (2-7r) ^ <?tt^ ( 1 1 ^ 1 1 2 ) = on R m . Hence, |x| ||) = on R m , which contradicts the condition 
ip(0) > 0. Thus, the function f m strictly decreases on (0, +oo) and, hence, / m (s) > / m (+oo) = for all s > 0. The 
Theorem \E\ is proved. □ 

5 Some statements involving a versatile general covariance function 

Previous results can be generalized to the class of positive definite functions built in [10] and used for the purposes 
highlighted in Section 1. 

Lemma 6. Let tp G C ([0,+oo) n ), n G N, and / °° . . . J °° |</?(ui, ■ • ■ ,u„)\ Yik=i u< k'^ 1 ■ ■ -du n < +oo for some 
dk G N, k = 1, . . . , n. Let h k ,b k G C(Ek), hk strictly positive in their arguments for all k = 1, . . . , n. Then 

K{x u ...,x n ,t 1 ,...,t n ) := v ,fJN!i > ...,JNi > ) f[b k (t k ) e m dl x...xR d -x£ 1 x...xi? n ) 

if, and only if, 

n 

gdi,...,d n (siy/hi(t 1 ),...,B n y/h n (t n ) y ) ]]_ b k (t k )(h k (t k )) dk /2 G $(-El X ... X E n ) 

fe=l 

for every Sfc > 0, k = 1, . . . , n, where 

9di,. 



poo pOO " 

!„(si,...,S n ) := / ... / ti£) TT Wfc fc_1 id fc /2-l(SfcWfe) d«l...<fti n - 

Jo Jo fe=1 
Proof. The statement can be proved in a similar way as Lemma [2j □ 

Let P_p be the set all finite nonnegative Borel measures on [0, +oo) n , n G N, and 

{/>oo />oo 
¥>(«!,...,««) =y ...jf e-(»«+^)# 1) ..,«J 1 , 1 6P» 

Obviously, = C[ , +O o) fl ^(o,+oo) and n£=i VfcC"*) 6 £ n f° r ever Y Vk G £\ k = 1, . . . , n. 
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Theorem 7. Let n G N. For all k = l,...,n, let E k be linear spaces, h k strictly positive functions such that 
hk G C(Ek) and d k G N. Then, the following statements are equivalent: 

1. K( Xl ,...,x n , ii,...,t„) := ^(jjf^, nLi {h k {t k )Y dk/2 G <f(R^ x ... xM d - x £ x x ... x E n ) 

^ e L". 

2. e -AMt*) e VA > 0, fc = 1, . . . ,n. 

Proof. Let us prove the implication (1) => (2). For every fixed k = 1, . . . , n in condition 1, we take </?(ui, . . . , u n ) = 
fkW), e £ x , and t, = £ S 4 for i ^ fc. Then Vk (jjf$|f) (M*fc))~ dfe/2 G $(K dt X E fc ) V^ fc G £ x . By 
TheoremHwe get e-^*^ G $(-E fc ) VA > 0. 

Let us now prove the reverse implication. Let e~ A/lfc ^ fc " ) G $(£'fe), VA > 0, fc = 1, . . . , n. By TheoremJU we have 
that <p k (jjf^f) (h k (t k ))- d */ 2 G $(R dfc x B fc ) V(^ fc G £\ k = l,...,n. We take (p k (u k ) = e-"*"*, w fe > 0. From 
definition of class £™ follows, that 9? (|f^, ■ ■ ■ , lff§) IlLi (M*k)) _<ifc/2 G $(R dl x . . . x R d - x E x X . . . X £ n ) 
Vv? G £ n . □ 
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